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We present a framework where first principles calculations of jet modification may be carried 
out in a non-perturbative thermal environment. As an example of this approach, we compute the 
leading order contribution to the transverse momentum broadening of a high energy (near on-shell) 
quark in a thermal medium. This involves a factorization of a non-perturbative operator product 
from the perturbative process of scattering of the quark. An operator product expansion of the non- 
perturbative operator product is carried out and related via dispersion relations to the expectation 
of local operators. These local operators are then evaluated in quenched SU (2) lattice gauge theory. 

PACS numbers: 



I. INTRODUCTION 

As of this time, the Large Hadron Collider (LHC) 
has completed 2 successful runs and there is a wealth 
of data on the modification of hard jets from the Rel- 
ativistic Heavy-Ion Collider (RHIC) [U E] and the 
LHC [i^HSj. With the similarity between the various 
soft observables between RHIC and LHC the study 
of jets has moved to the forefront of heavy-ion pro- 
grams at both these colliders. 

In the last several years, the science of jet quench- 
ing has undergone considerable evolution. There 
are now four different successful jet quenching for- 
malisms based on perturbative QCD (pQCD) [6]- 
[19] and a collection of formalisms based on 
AdS/CFT gUHlS]- While one would have ex- 
pected a large disparity between the physical pic- 
tures underlying the strong and weak coupling ap- 
proaches, there are actually considerable differences 
between the various pQCD based approaches. Be- 
sides the differences in the description of the per- 
turbative gluon emission process, the description of 
the medium is quite different in the various ap- 
proaches: In both the Armento-Salgado-Weidemann 
(ASW) and the Higher- Twist (HT) approach, one 
assumes that the transverse momentum exchanged 
in numerous interactions with the medium is soft 
enough that one may approximate the distribution 
as a Gaussian, and retain only the leading two mo- 
ments (mean and variance). The variance of this 
Gaussian transverse momentum distribution is of- 
ten referred to as q. In the Gyulassy-Levai-Vitev 
(GLV) formalism, the exchanged momentum is as- 
sumed to have a considerable hard tail, such that 
it cannot be approximated as a Gaussian broaden- 
ing. In the Arnold-Moore- Yaffe (AMY) formalism 
one describes the medium using Hard-Thermal-Loop 
improved perturbation theory [2S] . 

With the exception of the AMY formalism, none 
of the other pQCD based formalisms can be said 



to be a first principles calculation. In all cases the 
transport parameter q (either averaged or a normal- 
ized function of space-time in a fluid dynamical sim- 
ulation) is a fit parameter in the calculation, set by 
comparison to one data point. Even in the AMY for- 
malism the strong coupling constant as is varied to 
fit one data point. Thus, even the AMY formalism is 
not, strictly speaking, a first principles calculation. 
The strong coupling approaches, though first prin- 
ciples calculations, are not sufficiently sophisticated 
to address the great variety of jet modification data. 
The predictions from such calculations also seem to 
be inconsistent with the rising Raa observed at the 
LHC [26,. 

The goal of the present paper is to suggest a setup 
where a first principles calculation of jet modifica- 
tion can be carried out using a combination of per- 
turbative and non-perturbative methods. The per- 
turbative sector will be similar in form with the 
higher-twist approach in that it will involve a fac- 
torization of the perturbative sector describing the 
propagation of hard partons from operator products 
which will be used to describe the medium. The 
computation of these operator products in the non- 
perturbative sector will be carried out using finite 
temperature lattice gauge theory. We would point 
out already at this stage that a completely first prin- 
ciples calculation can never be directly compared 
with data. It will, however, provide constraints on 
the number, structure and normalization of the var- 
ious transport coefficients that one routinely uses to 
construct a phenomenological analysis of the data. 

The paper is organized as follows: In Sect. II, we 
describe the set up where calculations can be car- 
ried out and in particular we will attempt to justify 
why the current method to identify and estimate jet 
transport coefficients is the better alternative. In 
Sect. HI we will focus on the particular process of a 
hard quark propagating through a medium and set 
up the formalism for this process. In Sect. IV the 
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various regions of phase space will be explored. In 
Sect. V, dispersion relations that will be used to 
evaluate the operator products will be set up. In 
Sect. VI we discuss the details of the lattice gauge 
theory calculation. We conclude in Sect. VII with 
an outlook for future work. 



II. PQCD PROCESSES IN A QGP BRICK 

The notion that jet transport coefHcients repre- 
sent properties of the medium and thus should be 
calculable in lattice QCD has definitely been infor- 
mally considered for some time now. The most naive 
approach would be to simply take the expression for 
a given transport coefficient, say g, as derived in an 
appropriate effective theory in Ref [57], where 
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and attempt to compute this on the lattice (In the 
equation above Fj_'^'^ is a gauge field strength opera- 
tor, one of whose indices are either 1 or 2). This par- 
ticular form of the transport coefficient is obtained 
in either covariant gauge or light-cone gauge. 

The equation above is not manifestly gauge invari- 
ant and requires the introduction of Wilson lines. 
At first sight, the path taken by the Wilson lines 
seems arbitrary. However, following the arguments 
in Ref. [28], one obtains four different Wilson lines 
that need to be included, two along the light-cone di- 
rection and two along the transverse direction. The 
fully gauge invariant expression for q is now given 
as. 
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In the equation above, U represents a Wilson line 
along the (— ) light-cone direction and T represents 
a Wilson line along the transverse light-cone direc- 
tion. If the calculation were being carried out in co- 
variant gauge, only the light-cone Wilson lines will 
contribute, while for the calculation in light-cone 
gauge, only the transverse Wilson lines will con- 
tribute. Thus while the exact expressions are rather 
different in the two gauges, both may be derived 
from Eq. Given the extent of the Wilson lines 
(and the issues related with analytically continuing 



an euclidean operator product to one that is almost 
light-like separated) , it appears almost impossible to 
evaluate these on a finite size lattice. 

However, there exists an alternative, based on the 
similarity between q and the gluon distribution func- 
tion and the method by which parton distribution 
functions (PDFs) are evaluated on the lattice [291 - 
I32| . i.e., using the method of operator product ex- 
pansions. Imagine a high energy process e.g. the 
deep inelastic scattering (DIS) of an electron with 
momentum k off a single quark prepared with mo- 
mentum p, at one edge of a finite volume V which 
is maintained at a fixed temperature T ~ ^qcd- 
At this temperature the volume will be filed with 
strongly interacting matter, which at temperatures 
somewhat below Aqcd will be a hadronic gas and at 
very high temperatures will be quark gluon plasma. 
We maintain the chemical potential ^ — so that 
the contents have the conserved charges of the vac- 
uum. On scattering off the electron, the quark will 
produce a hard virtual quark which will then prop- 
agate through the medium. In vacuum such a par- 
ton would undergo a perturbative shower, spraying 
partons with ever lower virtuality until the scale be- 
comes comparable to Aqcd and hadronization be- 
gins to set in. In the presence of a strongly inter- 
acting medium the produced shower will scatter off 
the constituents in the medium, diffuse in transverse 
and longitudinal momentum, and be induced to ra- 
diate more partons leading to a further degradation 
in the energy of the part of the jet which escapes the 
medium. 

If the medium is not larger than E/^q, where E 
is the energy of the jet, and /io is the minimum scale 
below which pQCD is no longer applicable, a portion 
of the jet will hadronize outside the medium. The 
differential cross section for any particular outcome 
from such a hard scattering process can be expressed 
using the standard factorized formula, 



dah 



a 



k-pQ- 
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where C'^'^ is the usual leptonic tensor and dW^" 
is the differential hadronic tensor for the particular 
process of interest; all interactions which involve the 
QCD coupling g are contained within the hadronic 
tensor. 

Say further that in the hadronic tensor, we could 
factorize the initial distribution of the hard quark, 
the hard scattering off the photon and the final prop- 
agation through the medium as. 



(4) 



dW" = / dxf{x)da^"'D{{pf}). 



In the equation above, f{x) represents the distribu- 
tion of the initial quark; in the case of quark inside a 
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proton this would simply be the parton distribution 
function. In the case of a single quark it is simply 
(5(1 — x). The term a'^'^ represents the hard cross 
section for the scattering of a quark off a virtual 
photon. The function D which is a function of the 
set of measured final state momenta {p/} includes 
all final state effects after the hard collision of the 
quark with the photon. The general structure of D 
may be written as 

DiiPf}) = E(^l'^il^) 

X m]\{Pf}X){{pf}X\Q,\0), (5) 

where, \M) represents the medium where the jet 
interacts, |{pj}X) represents an inclusive hadronic 
state containing the detected hard momenta {pf} 
and other states which are not part of the medium. 
The operators Q, represent the part of the pro- 
cess which occurs outside the medium and fragments 
to yield the detected "non-medium" final state. The 
remaining operator Oj represents the part of the 
process that occurs within the medium. In a real 
heavy-ion collision such a distinction may be impos- 
sible to even formulate. However, in the theoretical 
scenario of a hard jet propagating through a finite 
medium, such a separation can be carried out order 
by order. 

In the case of a single inclusive measured hadronic 
momentum, D would become the standard fragmen- 
tation function (if Oj = 1 there would be no medium 
effect, otherwise one would obtain the medium modi- 
fied fragmentation function) . For more exclusive ob- 
servables (with more specified momenta), D would 
represent a more complicated object 33J. We should 
make it clear that the momenta which specify D do 
not need to be hadronic and may be completely par- 
tonic; in fact the particular D that we will consider 
will be completely partonic. 

In the remainder of this paper, we will consider 
evaluating Ot by perturbing in the weak coupling of 
the hard produced quark with the medium. Note 
that this does not assume that the coupling within 
the medium is perturbatively weak. We will encode 
the effect of the medium on the hard quark in terms 
of an infinite series of local, power suppressed, oper- 
ators (suppressed by powers of the hard scale Q^)- 
Thus Oi will be obtained as a series of local op- 
erators and ever more suppressed perturbative 
coefficients c!^/ [Q^] , 

= E J^^n- (6) 

While perturbation theory is valid for the inter- 
actions of the hard quark, it is not valid for the 



local operator products. Any evaluation in pertur- 
bation theory necessarily requires the specification 
of a gauge and the calculations in this paper will 
be no different. Each choice of gauge will result in 
a slightly different set of perturbative terms along 
with a slightly different set of local operator prod- 
ucts. For gauge invariant observables such as D the 
total sum will be gauge invariant. To demonstrate 
this however, one needs to be able to evaluate the 
operator products (for at least the first couple of 
terms) . 

In all prior attempts to evaluate D, the non- 
perturbative sector has never been evaluated exactly. 
In the HT scheme, which is closest in spirit to the 
present discussion, the operator products (or some 
combination of them) are treated as parameters of 
the theory. A model is assumed for how they would 
depend on intrinsic properties of the medium such 
as the temperature T. The overall normalization is 
set by comparing with one data point. In this paper 
we present the first effort to estimate these operator 
products non-perturbatively on the lattice. 

The primary motivation for this effort is to test 
if such an approach is at all feasible. There is no 
attempt to be exhaustive and only the simplest pro- 
cess of jet broadening will be considered: the broad- 
ening of a single quark by a single scattering with 
momentum exchange k± in a hot medium. Dividing 
the mean by the length of the medium will yield 
the transport coefficient q. The question that will 
be addressed in this paper is if such an approach 
is at all possible. To this end, we will calculate the 
perturbative part only in — gauge and the non- 
perturbative part in quenched SU{2) lattice gauge 
theory. In this sense, this paper should be viewed as 
a "proof of principle" of such a methodology. Issues 
related to renormalization on both the perturbative 
and the non-perturbative side will be ignored. The 
evaluation of the perturbative coefficient functions 
in an alternate gauge, the computation of the modi- 
fication of the shower pattern of the jet, and the eval- 
uation of the non-perturbative operator products in 
SU{3) will be left for future efforts. 

We note in passing that, while in this paper, we 
assumed the factorization of the hard scatting from 
the final state scattering, this (assumption) is not 
strictly necessary in such a framework. Indeed one 
may consider e+e~ annihilation within such an enclo- 
sure and calculate the modification of the back-to- 
back pair of jets. Depending on the choice of observ- 
able and gauge this will lead to a unique expansion 
in the form of Eq. ^ . 



III. LEADING ORDER DERIVATION 

In this section, the operator expectation D will be 
factorized into a pcrturbativc and non-pcrturbativc 
part. As pointed out above, we will consider the 
simplest process of jet broadening at leading order 
in the medium. To this end, we consider the propa- 
gation of a hard virtual quark through a hot medium 
with the quantum numbers of the vacuum. The large 
scale associated with this parton allow for the use of 
perturbation theory and we compute the first pcr- 
turbativc contribution which occurs only in the pres- 
ence of a medium. 

Imagine a quark in a well defined momentum state 
\q) = (j^, 0_l) impinging on a medium |M) and 
then exiting in the state 



\q + k) = 



[2{q-+k-)] 



,q +k ,k± 



with the medium state absorbing this change in 
momentum and becoming \X). The quark is as- 
sumed to be space-like off-shell with virtuality = 
29^(7" < with the negative z-axis defined as the 
direction of the propagating quark. In a physical sit- 
uation, one would have a gluon radiated off a quark, 
with either the gluon or the quark space-like off-shell 
(or both). The space-like parton would be placed 
closer to its mass shcill by scattering in the medium. 
The rate of scattering is controlled by the transport 
coefficient q. To mimic this process we have consid- 
ered the very simple process of a space-like quark 
scattering off the glue field in an extended medium. 
The case of an on-shell quark is included in the limit 
of Q2 ^ 0. 

Consider the reaction in the rest frame of the 
medium. In this frame f/' > 0, and we have defined 
the z-axis such that qz < 0. In this choice of frame, 
for a space- like quark we have q^ = {q^+qz )/V2 < 
and q~ = (g° — qz)/V^ > 0. If the z-axis were 
chosen such that > the g+ and q~ will sim- 
ply switch roles. For a space-like quark we have 
9" ^ \Qz\, and this implies that q~ > q'^ . For 
a jet one r equires q~ ^ g+. Alternatively stated 
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The spin color averaged transition probability (or 
matrix element) for this process, in the interaction 
picture, is given as 

W{k) = ^^{q-;M\T*e'fo ^tHAt) +k^^x) 

X {q- + k± , X\Te-' ^0 dtH,{t) \q-^M), (7) 

where, we have averaged over the initial color and 
spin of the quark, assuming that the medium is in 
a fixed state. In the case of a thermal medium 




FIG. 1: A quark scattering off a gluon in medium \M). 



one may use the density matrix to average out the 
initial state. We will assume that all this is im- 
plicitly included in \M) . In the equation above. 
Hi = J d^xi){x)igt''-i>'A1^{x)iP{x) and T (T*) rep- 
resents time(anti-time)-ordering. Expanding the ex- 
ponential to leading order yields, 



W{k) = 



2N, 



■{q-;M\ J dt 



xd^yi^iy) 4(y)V'(2/) 



X \q- +k±;X){q- + kj_;X\ 
X ilj{x)4{x)ilj{x)\q-;M), 



(8) 



where, = t^A^. To deal with the factors of time 
t and volume V, we introduce box normalization for 
the quark wave-functions and later take the limit 
of t,V — >■ 00. In box normalization, 'ip{x)\q~) = 
e~'''^u(g)/\/F, we get. 



W{k) = ^ I d'xd'yTr 



{M\^^ A{y) (9) 
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Shifting, the x and y integrations, the four volume 
may be extracted (/ d'^x = tV) and divided out by 
the factors in the denominator. The mean fc^ which 
yields q has the obvious definition. 



(10) 



where, we have summed over all values of the four 
vector k with the restriction that the final out going 
quark remain on shell. Where t represents the time 
spent by the hard quark in the thermal volume V. 
With the overall factor of four-volume removed we 
can take t,V ^ 00. 

We will now demonstrate that in the limit that 
q goes near on-shell, i.e., q~ ^ q~^, the expression 
above reduces to the well known expression for the 
transport coefficient q. Taking the limit that = 
2q~^q'~ — )• while q~ 00, we can simplify the 
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Dirac trace as 



{M\TryAU+ « A]\M) 
= 8{q-rTrrt>^]{M\At{y)A+{^)\M). (11) 

The imaginary part of the propagator yields the on- 
shell (5-function, which may also be simplified as, 
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Since, k~ has been ignored, compared to it may 
be integrated over to yield 27r5(y+). The fc^ may 
be combined with the vector potentials to yield, 
V_L^+ ~ F^. Absorbing both factors of k±, we 
obtain an expression containing only field strength 
tensors. 

Substituting the above simplifications, one ob- 
tains, 



X {M\F+-^iy-,y^)F+iO)\M). 



(13) 



This is the standard definition of q. Note that noth- 
ing is specified about \M), it may indeed be an ar- 
bitrary medium. If \M) is a thermal medium, then 
it must be averaged over in the sum over all ini- 
tial states. Averaging with a Boltzmann weight will 
yield. 
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Note that in the above derivation, no ordering is 
introduced between the two field strength operators. 
The expression above is not gauge invariant, but is 
gauge covariant. This implies, that if one were to 
carry out an operator product expansion in terms of 
local operators, one could reorganize the expansion 
to only contain gauge invariant local operators. Any 
gauge dependence would then only be contained in 
the coefficient functions. 



IV. THE OFF-SHELL REGIME AND THE 
NON-PHYSICAL REGIME. 

In the preceding section, we considered the pro- 
cess of a near on-shell quark propagating through a 
hot medium, at leading order in the scattering off 
the medium. In this section, the case of a slightly 
off-shell quark will be considered. The quark virtu- 
ality or offshellness will still be small compared to 
the energy. Once the operator products have been 



isolated, we will consider the process in the region 
of very high virtuality, of the order of the energy, 
and consider an expansion in a power series with 
increasing negative powers of the virtuality. 

Consider the imaginary part of the propagator in 
Eq. ([9]). In the limit that q^ is very large, and q^ is 
vanishingly small, there is a pole at the point where 
k~^ = {k'j_)/ {2q^). In the regime where q~^ ^ q~ 
but q^ is not vanishingly small (i.e., the parton has 
a non-negligible virtuality) we will obtain small ad- 
ditive contributions to the gauge covariant structure 
derived above. In this section we consider the more 
physical limit where q^q~ ~ fc^^ ~ X'^{q~)^, where A 
is a small dimensionless constant. In this case, the 
Dirac matrix structure will be simplified by taking 
the trace as, 

Tr[^ A{0)U+ « A{y)] - 4A^(0)G^.A''(y), 

with, G^, = [g^((7 + ky + q'^iq + kY - {q + k) ■ qg^,] 
Expanding this out, we obtain, 

A(0) • G ■ A{y) 

= 2q-A+{0)q-A+{y) + q-A+iO){q+ + k+)A-{y) 

+ q+A-iO)q-A+ + {q+ + fc+)A- (O)?" A+(y) 

+ q-A+{0)q+A-{y) + 2q+ {q+ + fc+) A" (0)^" (j/) 

- q-A+{0)kj_ ■ A^{y) - k^ ■ A^{Q)q-A+{y) 

- q+A-{0)k^ ■ A Ay) - ki_ ■ AAO)q+A-iy) 

- [q-iq+ + k+)+q+q-] (15) 
X [A+iO)A-iy) + A-iO)A+iy) - A^O) ■ AAy)]- 

We now consider this expression in A^ = gauge, 
where we may drop terms which go as Q'^/q~ ~ 
X^q~. This leads to a considerable simplification of 
the final expression, 

A ■ G ■ A = 2q- A+ {0)q- A+ (y) 
+ q-A+{Q)k^,^ ■ A^y) + • A'i{0)q-A+{y) 
- [q-{k+ + q+)+q-q+][A^A0)-A1iy)]- (16) 
The exponential phase factor is. 



exp 



2q- 



q ]y + k±,f,y'i 



(17) 



where the general (_L)-4- vector implies A± = 
[0, 0, Aj^]. Using these relations, we may simplify, 

2{q-r{-k^A±,^)A+{0)A+{y)e^'^(y^ (18) 
= -2(g-)2v^yl+(0)Vi,,,A+(y)e^^(^). 
The next set of terms simplify as, 

e'^q-A+{0)k^^,A'Ay)kl 

= 2{q-ftV^,^A+ [q+-zd+]AAA 

= 2(g-)2 [W^,^A+d+AAx)+iV^,^A+q+AAy)] 



(19) 
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The first term in the bracket above, can be com- 



bined with Eq. ( 18 ) to produce the field strength 
tensor at location x. There is another term similar 
to the one above which can be combined to form the 
field strength tensor at the origin. The last line in 
Eq. ([l6|) may be re-expressed as, 



^2{q- f [a+Ai,^(0)a+<(y) + 2zg+Ai,^9+A^(y) 
- ig+a+Ai,40)A^ + 2(q+)2Ai,^(0X(y)] .(20) 

The first set of terms in the equations above 
[Eqs.( |18|19|20| ] can be combined to obtain the 
known form that appears in the definition of the 
on-shell q, i.e. 2{q-fFl^^F'^^. Note that all 
terms in Eq. ( [20| are rather small [they scale as 
A22(g-)2v_L,;xA^0)V^y4+(y)] and thus the remain- 
ing terms may be ignored. 

We now have an expression for the transport coef- 
ficient q over a range of values of q^ where q'^ <C q~ 
but is still large enough that — 2q^q'' >• ^^cd- 
We can now take this particular operator product 
and consider its behavior over the entire complex 
plane of q'^ . 

We now analytically continue to the region where 
(7+ < and ^ ^ fc. Consider the analyti- 
cally continued, unphysical expression. 



Q 



Nc J (2^)4 72(7" 
(M|F+^(0)F+»|M) 
{q + ky+ie 



(21) 



We introduce a new object Q to indicate that the 
expression above is not the jet transport coefficient 
q. The discontinuity of the above expression in the 
region —q~ <C <C (?~ corresponds to q. 

In the regime where q^ ^ q^ ^ k, one can expand 
out the denominator as. 



1 

(Q2 - fc^ + 2g • fc) ^ Q2 2^ 



^^^1.(22) 



The instances of the gluon momentum k may be re- 
placed with derivatives. Adding, gluon scattering 
terms, we can convert the regular derivatives into 
covariant derivatives. Thus we obtain a series of 
gauge covariant expressions for the jet transport co- 
efficient. 



^ _ An^a, f d^yd^k ,,.yV2q- 
~ - ' (27r)4 ^ Q2 



{M\F+^iO)J2 



(23) 



n=0 



-q-iV-VV, , , > 



With all instances of k removed from the integrand 
(except for the phase factor), the integrals over all 



components of k can be carried out to yield four 
(5-functions over the position y. This yields a very 
simple expression for q in = gauge, in terms of 
local gauge invariant operators, 



4V27r2 



asq 



{M\F+' 



E 

n=0 



F+ IM).(24) 



The above expression requires some discussion. 
The discontinuity in the expression above across the 
real axis of corresponds to the transport coef- 
ficient q when —q^ <^ q~^ q~ . For g"*" ~ q~ 
and positive, there is another source of a disconti- 
nuity, from real hard gluon emission. This part is 
perturbatively calculable as long as = 2q~^q~ 
^QCD ^'^d does not depend on any properties of the 
medium. In the region where is space-like or 
g+ ^ ~^QCD there is no discontinuity across the 
real axis. Alternatively speaking, in the deep space- 
like region the internal quark-line cannot go on-shell. 
For virtualities which are not in the deep space-like 
region, the quark can still absorb a gluon from the 
medium and go on-shell and there will be a discon- 
tinuity. 



V. DISPERSION RELATIONS 

In the preceding section, the expression for q was 
generalized to the region of (a physically realizable) 
non-zero virtuality and then considered in the re- 
gion of (unphysical) very high virtuality. In the cur- 
rent section the two expressions will be related via 
dispersion relations in the complex q~^ plain. The 
expansion in the unphysical region [Eq. ( 24 )] will be 
used to estimate the value of q in the physical region. 



for 



X^q- 



In order to evaluate q — Disc Q 

we will use the method of dispersion relations: We 
will evaluate a similar integral in a region of the q^ 
complex plain where there is no discontinuity and 
use methods of contour integration to relate the eval- 
uated integral to q. 
Consider the integral, 



dq^ Q{q+) 



27ri {q^ 



(25) 



where Qo is large and positive. The contour is taken 
as a small counter-clockwise circle around the point 
-Qo- The residue of this integral is given as. 



1+ — 



:Q(g+) 



(26) 
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While this analysis can be carried out for arbitrary 
TO, we consider, for definiteness, the case of to = 1. 



In the limit where ^ XQ, we obtain Eq. (24 1 
with replaced by —2q^Qo, i.e. 



/i 



71 — 0^ ^ 



.(27) 



Since ,q 
simplified as. 



3> fc^, the above operator relation in 



Na2Qo 

Nc2Qo 
-iV 



n=0 



2Q' 



-w- 



n— m— 



2g' 



2r 

4%/27r2Q;, 
iV,2Qo 

(TO + fc)! f-iV 



\M). 



OO 

E 

A:=0 



fc! 



m=0 



2Q0 



2r 



(28) 



We can now deform the contour and evaluate it 
over the branch cut from q'^ > — A^Q to g+ — > oo. 
This yields. 



h = 



dq+^L]^e^^y- 



S k+ 



2q- 



Nc J " (27r)4 
{M\F+^^{0)F+{y)\M) 



2q' 



(9+ 
dq^ 



-Qo) 



Qo 



+ / dq+Viq+). (29) 



The second term in the equation above, refers to 
the contribution to the operator above from vacuum 
gluon radiation, i.e., the Bremsstrahlung radiation 
of gluons from an off-shell quark. As such, it con- 
tributes only in the region where the virtuality of 
the incoming quark is time-like and is independent 
of the temperature of the medium. Thus for a fixed 
T the second term above is a constant, while the first 
depends on the temperature of the medium. 

The limits on the first integral in the equation 
above allow for a simple expansion of the denomi- 
nator. The factor Q*^ ^ Q is much larger than the 
q^ ^ X^Q in this region and thus we obtain, the 
much simplified relation, 



dq' 



Qo 



OO 

E 

n=0 



Qo 



h 



dq+V{q+)m 



To obtain q, a general functional form in the vicin- 
ity of —X^Q < q^ < X^Q must be used. We start 
with the assumption that g at a fixed is a slowly 
varying function of q^ . This allows us to use a trun- 
cated Taylor expansion for q [We should point out 
that using the first few terms of the Taylor expansion 
is, in itself, an assumption regarding the functional 
form of g(g''')]. To provide a simple illustration of 
the procedure, we take only 3 terms; in the final nu- 
merical results we will only use those results where 
the first term greatly dominates over all subsequent 
terms (Note that an arbitrary number of terms in 
the Taylor expansion may be retained for a more 
accurate determination of q), 



Qiq ) = q + qq + 



q"{q+Y 



(31) 



In the above equation q' = dq/ dq^\q+^Q. 

Using the above truncated Taylor expansion we 
obtain. 



h 



Q+ 

Jdq-^ 

-Q+ 



q QO ^ q 2 



+ / dq+V{q+) = 



g'2(Q+ 
SQl 





+ ^3 



. 2(Q+) 
' 3Qf, 



q"iQ+r 
3Qo 

' 5QI ■ 



(32) 



In the equation above, Q'^ represents the limit of 
integration over q^ for the jet. For a jet with max- 
imum virtuality fi^ and (— ) momentum q~ , Q^ = 
/{2q^). One may now simply compare with the 



expression for Ii from Eq. ( 28 1 and equate the vac- 



uum subtracted coefficients of Qq . 

The methodology outlined above can be made 
even more precise and straightforward by setting a 
definite value for = q~ . While this will readjust 
the relative importance of the various terms in the 
series it allows for simpler set of operators that need 
to be evaluated numerically. This simplifies Ii in 
Eq. pSl) to. 



and similarly simplifies Eq. (32) with Qq replaced 
by q^ . For a virtuality fi^ such that Aq^^ <C /^^ <C 
((7^)^, we can define a q^ or virtuality averaged q 

as. 



q{Q+)2Q^ 



-Q+ 
2qQ+ 



dq+q{q+) 
q"{Q+f 



(34) 
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where the second hne is only vahd in the hmit that 
g is a slow function of (or alternatively stated 
(5+ q^). We can obtain an estimate of this by 
studying the 2nd term in Eq. ( 33 1 . If this term is 



comparable to the first term then the above approx- 
imation is no longer valid. If this term is small, then 
one may obtain a good estimate of q from just the 
first term in the series in Eq. ( 33 1 . In the subse- 



quent section the forms of the operators and their 
evaluation on the lattice will be discussed. 



VI. LATTICE CALCULATIONS 

In the preceding sections, the jet transport pa- 
rameter q, as obtained in the physical regime of jet 
momenta q^ ^ X^q^ ^ q^ , was related via disper- 
sion relations to a series of local operators in an un- 
physical regime where g"*" = —q~. The availability 
of a series of local operators, suppressed by powers 
of the hard scale q~ allow for the calculation of such 
non-perturbative operator products on the lattice. 
In essence, our task is to compute the finite temper- 
ature Minkowski space correlator. 



?i|e-^^Oi(i)O2(0)|n) 



(35) 



in the limit that t 0. In the equation above, /3 is 
the inverse temperature (/? = 1/T), H is the Hamil- 
tonian operator, and \n) represents an eigenstate of 
the Hamiltonian. Using the standard relations of the 
imaginary time formalism of finite temperature field 
theory, we can relate the Minkowski correlator with 
the Matsubara correlator in Euclidean space. 



-IT) 



A(r) = TV 



-JdrHir) 



O,(r)O2(0) 



,(36) 



for the case where there are no time derivatives in 
Oi and ©2 and yields, V>{-iT) = i^*A(T) for a 
total of Nt time derivatives in D^{t). As a result we 
obtain the simple relation that. 



2?>(i = 0) =^^*A(T = 0). 



(37) 



Using the above relation, the local operator products 
in Minkowski space may be obtained from the local 
operators in Euclidean space. 

In the following, we list out the operators that 
must be evaluated and re-express them in a form 
where they may be easily calculated on the lattice. 
In this first exploratory attempt, the calculation will 
be carried out for an SU (2) gauge theory on a space- 
temperature lattice in the simplified quenched ap- 
proximation. Quark-less SU{2) possesses a nega- 
tive /3-function as in full QCD. Since issues of 



higher order contributions and renormalization were 
ignored in the perturbative sector, renormalization 
will be ignored in the non-perturbative sector as 
well. The extension to more sophisticated simula- 
tions in quenched (or unquenched) SU{3) will be left 
for future efforts. In defense of the current effort, we 
point out that in the context of jet transport coeffi- 
cients in heavy-ion collisions, quenched calculations 
may provide a very realistic estimate, as the early 
dense plasma is believed to be gluon dominated. 

In the language of links, the field strength tensor 
f^F^^^ may be expressed as. 



F., 



U - W 

2igaL 



(38) 



Where, C/^^ represents a plaquette in the fii' plane 
and fli is the lattice spacing. Similarly, terms with 
a covariant derivative may be expressed as. 



-DiF^^ix) 



Ff,t,{x^+aL,x) - U4{x'\x)Ff^^{x*x) 



ah 



,(39) 



where, 6^4 represents a gauge link in the 4-direction. 
In this paper, we have only used the right derivative 
as we seek only an order of magnitude estimate of 
terms with a time derivative, as argued below. 
The first operator to be evaluated is 

{M\Fl^F^ + \M) = Y,{n\e-f'"Fl^F^'+\n)m 



E 



e-f^"{n\Fl^F^-+\n) 



where, |n) represents an eigenstate of the full Hamil- 
tonian. We do not indicate the location of the two 
F field strength tensor insertions as both are at the 
same location. 

We now discuss the rotation of the operator prod- 
ucts to Euclidean space. This involves the two rota- 
tions: 



-ix^ and A° ^ iA^ 



F" 



iF^ 



(41) 



As a result, 

^ (F^^F^^) - {F^^F^^).{42) 

In the above equation, we have ignored terms which 
measure topological charge density (F^^F'*^), as- 
suming instanton effects to be small. 

We report results on a (4 x nt)"^ x rtt lattice where 
nt is varied from 2 to 6. Note that finite temperature 
calculations are meant to be carried out in the limit 
that nt <ti rig. For = 2 we have also carried out a 
calculation with ng = 12. We have not repeated the 
calculation with smaller values of nt and the largest 
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value of Ug = 24 as the results for nj = 2, 3 do not 
seem to be have any dependence on for > 12. 

For this first attempt we will use the Wilson gauge 
action for SU (2) ^34*, "351 . The scale (or lattice spac- 
ing) is set on the lattice using the renormalization 
group formula [511 [55] 



1 



V247r2 J 



exp 



12^ 



(43) 



where, is the lattice spacing, g represents the 
bare lattice coupling and represents the one 
dimension-full parameter on the lattice. Compar- 
ing with the vacuum string tension, we have used 
Al = 5.3 MeV. For a lattice at finite temperature or 
one with nt <^ ris, the temperature is obtained as 



T 



1 



(44) 



While this may not provide the best means to set the 
scale on the lattice at finite temperature, it proves 
sufficient for a "proof of principle" approach as out- 
lined in this paper. Gauge configurations are gener- 
ated using a simple heat bath algorithm [52]. Cal- 
culations with nt = 2, 3 consist of 1000 heat bath 
sweeps for each data point, those for rit = 4, 5 rep- 
resent 3000 sweeps and the calculations with rit — 6 
were done with 5000 heat bath sweeps. 



□^nj,=3, nj=12 
11^=4, nj=16 
A^iif =S, ns=20 



I ■ I ■ I 




100 200 300 400 500 600 700 800 



T (MeV) 



FIG. 2: The temperature dependence of the local op- 
erator {f^'F"*"') scaled by T* to make it dimension- 
less. Given the expression for the lattice spacing in 
Eq. ( |43[ l a phase transition is expected in the vicinity 
of T = 150 - 200 MeV. This is also reflected in the ex- 
pectation of the operator product. See text for details. 



ries [of Eq. (33)] as a function of the temperature as 
measured on the lattice. We find that while the cal- 
culations with nt = 2, 3 do not show scaling with lat- 
tice size, the calculations with rit = 4, 5, 6 show very 
good scaling especially in the region away form the 
phase transition, i.e., above T — 200 MeV. With the 
particular value of quenched (quark- less) SU{2) 
lattice gauge calculations demonstrate a phase tran- 
sition in the region between 150 — 200 MeV. This 
is also borne out by the expectation of the operator 
product. We scale the results with to construct 
a dimcnsionless ratio. 

For this first attempt to extract q from a lattice 
calculation we have not carried out a rigorous esti- 
mate of the statistical error. For the data points 
at T > 300 MeV, for > 3, the statistical er- 
ror is approximately the size of the plotting symbols 
(this may be easily surmised from the scatter in the 
curves). Estimates of the statistical error for oper- 
ators composed solely of gauge fields in quenched 
SU{2) have appeared in several papers e.g., see 
Ref. |36j . In our view, the more important error 
is the systematic error brought on by the neglect of 
the the higher derivative terms in Eqs. ( 32|33 1. To 
estimate the value of q solely from the first term in 
the expansion in Eq. ( 33 1 on requires that the higher 
derivative terms be small. As an estimate of the size 
of these terms we compare the modulus of the ex- 



pectation of the first two operators in Eq. (33), for 
the case of rif = 4 in Fig. |3] 



n=4, n =16, q" = 20 GeV 

t ' s ' ^ 




100 200 300 400 500 600 700 

T (MeV) 

FIG. 3: The temperature dependence of absolute values 
of the the local operator (F^'F^^) and the second oper- 
ator product {[F+^id'^F+'] /q~) . See text for details. 



In Fig. [2] we present our results for the calculation 
of the operator product in the first term of the se- 



The plot in Fig. |3] demonstrates that for temper- 
atures below T = 600 MeV, the expectation of the 
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operator /q~ for a g~ ^ 20 GeV is less 

than 25% of the first operator product. Thus for 
temperatures below 600 MeV, for a,q^ ^20 GeV, we 
may obtain an estimate of the transport coefficient q 
from only the leading term in this lattice calculation. 
As pointed out in the discussion of Fig. [2j the sta- 
tistical error is small enough only for T > 300 MeV 
in order to extract a continuum limit. Thus one can 
extract q from such a calculation only in the range 
300 MeV< T < 600 MeV. This range coincides with 
the highest temperatures reached at RHIC and LHC 
and thus will allow future, more sophisticated, ef- 
forts to compare meaningfully with the values of q 
obtained from phenomenological analysis of RHIC 
and LHC data. This constitutes the primary result 
of the current manuscript: the demonstration that 
the framework developed in Sections HI, IV and V 
can be used to obtain reliable estimates of jet trans- 
port coefficients in a hot medium. Of course, com- 
parisons with experiment will require both a more 
sophisticated perturbative analysis as well as a much 
more developed lattice calculation. 



VII. ESTIMATING q AND CONCLUDING 
DISCUSSIONS 



In this concluding section, we attempt a simple 
minded extraction of the jet transport coefhcient q 
from the lattice calculation outlined above. Recall 
that our calculation required that the hard quark 
moves through the medium without undergoing any 
radiation. This constrains the highest virtuality that 
the quark may possess for such an approach to make 
sense. In a future effort, partons with a higher ini- 
tial virtuality will be considered. These will undergo 
radiative splitting in the medium and may show sen- 
sitivity to a somewhat different series of operator 
products. 

We choose the 3rd last point in the 
curve in Fig. |3] This corresponds to a tempera- 
ture of T = 363 MeV, which is very close to the 
top temperature reached in RHIC collisions. At 
a T = 363 MeV, {F+'F+') = 0.04 GeV''. Also 
we are considering a lattice with a length given by 
4 X ntttL = 4/0.363 GeV^ = 11 GeV"!. This 
states that the maximum virtuality of a jet (with 
a. q~ — 20 GeV) which traverses such a length with- 
out undergoing radiation is given as — E/L — 
20/11/^ ~ 1.3 GeV'^. Thus Q+ = 1.3/40 GeV. 
With these estimates, we obtain. 



quenched SU{2) plasma. In most phenomenologi- 
cal estimates one quotes the q of the gluon. If the 
above calculation were done for an SU{2) gluon, 
the q would differ only by the overall Casimir fac- 
tor of Ca/Cf = 2N^/{Nl - 1) = 8/3 yielding a 
qc = 1.432 GeVVfm, at a T = 363 MeV. 

In future efforts, the calculation will be extended 
to higher statistics runs to evaluate the q across the 
phase transition. The next step is to evaluate the 
required operator products for a realistic jet which 
starts at a higher virtuality and undergoes radiative 
splitting in the medium. In such a calculation, the 
range of operators that will need to be evaluated 
[i.e., number of terms in the series of Eq. (33 1 that 



need to be retained] will depend on the particular 
parton in the shower, in particular on that parton's 
energy and virtuality. 

Finally, to be of use to jet modification at RHIC 
and LHC the calculation will have to be extended 
to unquenched SU{i), renormalization factors will 
have to be introduced, and better means to set 
the lattice spacing will have to be used. At this 
stage we may only set suggestive limits on such a 
future estimation: the quenched SU{2) calculation 
has 3 colors of gluons as the fundamental fields in 
its Lagrangian, whereas there are 8 colors of gluons 
in quenched SU{3), along with 3 colors of quarks 
and antiquarks in the unquenched SU (3) calculation 
(Note that even if the plasma were completely per- 
turbative, quarks would contribute differently to the 
calculation of q than gluons [38], or rather the lat- 
tice calculation could change considerably with the 
introduction of dynamical fermions). Ignoring such 
subtleties, and assuming a simple scaling law with 
number of fields in the Lagrangian, we estimate that 
the full q at RHIC would lie in the range: 

q{T = 363MeV) = 3.7GeVVfm - 6.5GeVVfm.(46) 

We should point out that while the above estimate 
is very specific to a particular range of , q^ of 
the propagating parton, the estimate obtained form 
phenomenological analysis of RHIC collisions is an 
average over a wide range of parton energies and vir- 
tualities. In spite of the many shortcomings of the 
above calculation, we find the very encouraging re- 
sult that our estimate for q at T = 363 MeV is com- 
parable with that extracted from phenomenological 
analysis of RHIC and LHC data [Ml HO]. 



(45) 



Using a^(1.3GeV^) = 0.375 [37], we obtain q = 
0.537GeV^/fm for an SU{2) quark traversing a 
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